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Abstract 

We show that Kompaneetz equation describing photon diffusion in an 
environment of an electron gas, when linearized around its equilibrium 
distribution, coincides with the relativistic diffusion discussed in recent 
publications. The model of the relativistic diffusion is related to sol- 
uble models of imaginary time quantum mechanics. We suggest some 
non-linear generalizations of the relativistic diffusion equation and their 
astrophysical applications (in particular to the Sunyaev-Zeldovich effect). 

1 Introduction 

An evolution in time of a stream of particles can in general be described by 
the Boltzmann equation. In this paper we are interested in an evolution of a 
relativistic gas of massless particles . The relativistic Boltzmann equation is 
the standard tool in a study of the evolution of the photon in a gas of charged 
particles [1][2]. A diffusion approximation to the non-relativistic Boltzmann 
equation has been applied in many models of the transport phenomena [3] . The 
diffusion approximation to the relativistic Boltzmann equation is applied in the 
heavy ion collisions and the theory of quark-gluon plasma [4] [5] . A long time ago 
Kompaneetz [6]( see also [7]) derived a non- linear diffusion equation applicable 
to the photon propagation in an environment of an electron gas. The equation 
describes an evolution of a beam of photons disturbed by Compton scattering 
and Bremsstrahlung. These are the processes widely studied in astrophysics 
[1][8][9]. 

There seems to be the unique way to approximate the non-relativistic trans- 
port phenomena by the diffusion. The concept of a relativistic diffusion did 
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not achieve such a generally accepted consensus. Various versions of the rela- 
tivists diffusion are discussed in [10] [11] [12] [13] [14] [15] [16] [17] (for a review and 
further references see [18] [19]). In refs.[10] [11] the relativistic Brownian mo- 
tion is uniquely defined by the requirement that this is the diffusion whose four 
momentum stays on the mass-shell and the generator of the diffusion is the 
second order relativistic invariant differential operator. Starting from this basic 
assumption we discussed in [20] friction terms leading to an equilibrium. In [21] 
we studied a diffusion of elementary particles with spin. 

In this paper we show (sec. 3) that if the Kompancctz equation is expanded 
around its equilibrium ( Bosc-Einstcin) solution and non-linear correction terms 
are neglected then the resulting diffusion equation is the same as the relativistic 
diffusion with friction . First, in sec. 2 we formulate the relativistic diffusion 
equation for masslcss particles in a form appropriate for a comparison with 
the Kompaneetz equation. In sec. 4 we show that the linearized Kompaneetz 
equation can be solved applying methods of quantum mechanics because the 
relativistic diffusion can equivalently be described as quantum mechanical evo- 
lution in an imaginary time. Finally, in sec. 5 we suggest a generalization of 
the Kompancctz equation to a moving frame and a background gravitational 
field. This is a non-linear equation taking into account the quantum statis- 
tics of bosons whose linear approximation coincides with the (linear) relativistic 
diffusion discussed in [17] (without friction) and in [22] (with friction). 



2 The relativistic diffusion of massless particles 

Following Schay[10] and Dudley[ll] we define the generator A m of the relativis- 
tic diffusion as the 0(3, 1) invariant second order differential operator on the 
mass- shell Tim 

2222222 \ 

P =Pa -Pi ~P2 -P3 = m c ■ (!) 

The limit m — > of A m is divergent. In order to obtain a finite result we take 
the limit m — > of m 2 A m . This limit coincides with another definition of the 
diffusion [21] as the process generated by M fiU M' 11 ' where M M „ are generators 
of the Lorentz group (the formula for a diffusion of a massless particle has been 
also derived in [23] as a continuum limit of discrete dynamics of space-time, see 
also [24]). 

Choosing p (boldface type letters denote three dimensional vectors) as co- 
ordinates on Ti.Q we obtain as the limit m — ► of m 2 A m the operator 

A H = P 3 Pkd'd k + 3 Pk d k (2) 

here k = 1, 2, 3 and d- 7 = The relativistic diffusion equation in the momen- 
tum space (a relativistic analog of the Brownian motion) reads 

d T <Pr = \l 2 A H c^ T . (3) 
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7 is a diffusion constant. The photon as a massless particle is a typically rel- 
ativistic object. At the same time this is a quantum system with an internal 
angular momentum (spin). In this paper we neglect the spin. In [21] we discuss 
the diffusion of particles with a spin (the helicity in the massless case). It is 
shown in [21] that if m = then the dissipative part of the evolution does not 
depend on the helicity . Hence, the neglect of spin in this paper is justified. 
We add a drag term 

^ ^ dpj + ^ ^ j dpj + P>1 dXfj, 

to the diffusion (3). The vector field R defines a perturbation (friction) of the 
relativistic dynamical system (at 7 = 0) 

^ = Y (4) 

dx " K ^ 
" ^ (5) 

where [i = 0, 1, 2, 3 and from eq.(l) (for m = 0) we have p° = |p|. r is an affinc 
time parameter describing relativistic invariant dynamics. Its relation to the 
laboratory time x° follows from the zeroth component of eq.(5). We set Y = 
till sec. 5 where the vector Y will describe a geodesic motion on a manifold. 
Let 

Q=^ 2 A H + X. (6) 
Now, the diffusion equation reads 

d T <t>T = G^t- (7) 

We are interested in the behavior of the relativistic dynamics at large values 
of t. In particular, whether the diffusion generated by Q can have an equilibrium 
limit. We say that the probability distribution c?xdp$ is the invariant measure 
for the diffusion process if the statistical mean values are r-independent (there 
may remain the dependence on x°) 

J dxdp<&(x, p)(f> T (p, x) = J dxdp$ T (p, x) </>(p, x) = const. (8) 

Here, the evolution of <& is defined by 

a T $ T - g*$ T (9) 

and Q* is the adjoint of Q in L 2 (dpdx) (the Hilbert space of square integrable 
functions;the Lebesgue measure dxdp is not Lorentz invariant, the factor | p» | 1 
necessary for the Lorentz invariance is contained in see [20] [22]). 
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The r-invariance (8) is equivalent to 

e*$ = o. (10) 

Eq.(10) defines a transport equation in the laboratory time x°. 

We denote by &er the x-independent solution of eq.(10). Then, R can be 
expressed in terms of <&er 

R i = \pj + \pjPkd k In <$>er- (11) 

We assume that <&er is a function of the energy p c multiplied by a constant 
of the dimension inverse to the dimension of the energy (/? = in the 
conventional notation, where T denotes the temperature). In such a case from 
eq.(ll) 

Ri = \pj + ^ C ft|p|(ln$ £ «)'( C /3|p|). (12) 

Now, eq.(7) reads 

d T (t> T = ^ 2 A H <p T + \~f 2 P] (l +pc\p\{ln<S> ER y)&cj> T . (13) 

Let us consider an initial probability distribution of the form $ = Qer.^- Then, 
from eq.(9) we obtain its evolution 

*r = ^er^t (14) 

as a solution of the equation 

(15) 

(derivatives over the space-time coordinates have an index x,till sec. 5 we re- 
strict ourselves to the diffusion in momentum space, hence we drop the spatial 
derivatives). The diffusion generator (2) is degenerate. This is easy to see if we 
express it in the spherical coordinates on Ho 

p =r (16) 

Pi = rcos(j)sm9, p 2 = rsin<psm8,p 3 = rcoa8. In these coordinates 

A* = r a &+3r& = l 6 r + 2& (17) 

where we introduced an exponential parametrization of r 

r = expit (18) 

Here, u varies over the whole real axis. It follows that the three dimensional 
diffusion of massive particles becomes one dimensional in the limit m — > 0. 
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In the coordinates (18) it is sufficient if we restrict ourselves to the drags 



X = ^Rl (19) 

The diffusion (7) is generated by 

g=\{d 2 u + 2d u ) + ^Rd u . (20) 

We may write for integrals of spherically symmetric functions 

dp = 47r<iwexp(3u). (21) 

Then 

g* = e- 3u g+e 3u , (22) 

where 

g + = ^{di-2d u )- 1 2 d u R (23) 

is the adjoint of g in L 2 (du). The probability distribution evolves according to 
eq.(9). The invariant measure dxdp&n = dxdu^i solves an analog of eq.(10) 

g*$> R = g + ®i = o (24) 

where 

= cxp(3w)$ K . (25) 



3 Linearization of the Kompaneetz equation 

In this section we compare the probability distribution (9) resulting from the 
relativistic diffusion equation with the density of photons moving through an 
electron gas and described by the Kompaneetz equation. We assume that &er 
in eqs.(14)-(15) is the Bose-Einstein equilibrium distribution. This is also the 
equilibrium solution of the Kompaneetz equation. We show that if the density 
of photons is low (so that the quantum bunching effects can be neglected) then 
the relativistic diffusion equation (9) and the Kompaneetz equation linearized 
around the equilibrium coincide. We begin with the relativistic diffusion (9). 
The x° independent solution of eq.(24) is denoted &ei- The drift R is related 
to <E>_ej (it is defined by $er in eq.(12)) 

R= i(3„ln$ B/ -2). (26) 

If the diffusion (9) starts from an initial distribution 

= (27) 
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then using the relation (26), we obtain the evolution of ^ 

d T ^ T = i 7 2 W + ^ 2 (d u ln^ EI )d u ^ T . (28) 

First, let us consider Jiittner equilibrium distribution [25] in L 2 (du) (the 
classical limit of the Bose-Einstein distribution) 

® J EI = r 3 exp(-/3cr) = r 3 n,j. (29) 

Then, from eq.(26) 

R=~(-l + 0cexp(u)) (30) 
The diffusion generator corresponding to the Jiittner distribution reads 

S = \l 2 (d 2 u + Zd u -(3cexv{u)d u ) 
= \ 1 2 (r 2 d 2 r +Ard r -f3cr 2 d r y 

For the Bose-Einstein distribution (in Compton scattering the photon number 
is preserved; if an equilibrium is achieved, then the chemical potential /i ^ 0, 
see [7]) 

$1, = r 3 ( exp(/3( M + cr)) - l) = r 3 n E (n) (32) 

we have 

R = ~7j7 2 ( — 1 + P c cxp(u)^l — exp(— /3/U — /3c expu)^ ^. (33) 
Eq.(28) for an evolution of (with the drift (33)) in the r-coordinates reads 
T tf J = i 7 2 r" 2 (9 r r 4 a r * 7 + r 4 (<9 r lnn s )<9 r * / ), (34) 

where 

d r lnra s = -/3c(l -exp(-/3 M -/3cr))- 1 . (35) 
In the high energy limit we have in the diffusion equation (34) 

d r lnn^ — ► —fie 

(the same as for the Jiittner distribution). 

Now, we can compare the relativistic diffusion equation with the (non-linear) 
Kompaneetz equation usually written in the form [6] 

d T n = K 2 p- 2 d p (p 4 (d p n + n + n 2 )y (36) 

In order to explore the appearance of friction in eq.(36) we expand the 
photon distribution of the Kompaneetz equation around its equilibrium value 
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in the same way as we expanded the distribution of the diffusion process (15) 
around the equilibrium in eq.(28). If we neglect n 2 on the rhs of eq.(36) then the 
rhs disappears (because d p nj = —nj, with p = (icr) for the Juttner distribution. 
Let 

n = exp(-p)x- 

Then, neglecting the n 2 term in the Kompaneetz equation (36) we obtain (after 
an elementary rescaling r — ► (icr = x, k 2 = ^j 2 ) the diffusion equation (34) 
with tie — > nj. 

In general, the Kompaneetz distribution n equilibrates to rig (32) (as d p riE + 
tie + n% = on the rhs of eq.(36)). In more detail, let us write the initial 
condition for the Kompaneetz equation in the form 

n = n E x 

then the evolution of \ is determined by the equation 

drXr = K 2 p- 2 (d p p i d p x T + P A {dp hi n E )d p Xr ) + K 2 p- 2 n^d p (p A n 2 E {jdr ~ Xr)) 

(37) 

The last term in eq.(37) is absent in eq.(34). It describes the bunching effect 
of the quantum statistics which promotes bosons to condense. As long as the 
density of photons is low this term is negligible. 

The Kompaneetz equation finds applications in astrophysics [1] [26] [27] [28] [29] [30] 
[31][32][33]. In particular, the CMBR photons described by tie{p — 0) are scat- 
tered when passing clusters of galaxies and other reservoir of hot plasma. When 
applied to the Sunyaev-Zeldovitch effect the Kompaneetz equation is usually 
discussed only in a linear approximation. In such a case 

tt T (r) = / dyy-ip T (r, y)*(y) = (27r 7 2 r)^ / °° f cxp ( - ^(ln f - f 7 2 r) 2 ) 

(38) 

is the solution of the equation 

a T * T = i 7 2 r- 2 a r r 4 a r i' r = g*t> T 

with the initial condition The solution (38) is discussed in [27] [32]. For small 
r and the initial condition ^ — ue{p — 0) the solution can be approximated by 

* r = nE(p = 0)+Tg*n E (p = 0). 

The discussion of Sunyaev-Zeldovich CMBR spectrum distortion [27] [30] [31] [32] [34] 
is usually restricted to this approximation. The Kompaneetz equation has been 
derived from the Boltzmann equation [6][1] under an assumption that the elec- 
tron velocities are non-relativistic. The diffusion limit of the Boltzmann equa- 
tion with relativistic corrections to the Compton scattering have been calculated 
in [9] [31] [35]. 
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4 A relation to the imaginary time quantum me- 
chanics 

In this section we show that the (linear) relativistic diffusion equation can be 
solved by means of the methods of quantum mechanics. We can express solu- 
tions of the diffusion equations by an imaginary time evolution generated by a 
quantum mechanical Hamiltonian.The diffusion generator is of the form 

G = \l 2 dl - Lud u (39) 

( we denote a general drift by u here in order to distinguish it from particular 
R of the earlier sections) . Let 



Q(u) = / lu. (40) 



Then 



where 



exp(-0)g cxp(ft) = -H = ]p 2 dl - V (41) 
V = \co 2 - \d u w. (42) 



Let 



ip = exp(-fi)<£. (43) 

Then, from eq.(41) it follows 

4> T (u) = (exp(rg)0)(u) = exp(fi(u))(exp(-rfl")^)(u). (44) 

Inserting (f> = 1 in eq.(44) we obtain exp(— tH) exp(— 0) = exp(— O). Hence, 

Hexp(-Q) = 0. (45) 
For the Jiittner model when 



eq.(42) gives 



V = h 4 + l 1 ^ 2 ^ 6Xp(2u) ~ I (1 + 3 7 2 )7 2 /3cexp( u ). (46) 

This is a soluble model (the Morse potential ). The eigenvalue equation for 
H has a solution in terms of hypergeometric functions [36]. The propagation 
kernel determining the time evolution (44) is also known. The eigenstate of H 
with the zero eigenvalue is 

exp(-fi) = cxp (Ju - i/?ccxp(u)) . (47) 
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The probability density in the state (47) 

exp(-2ft(u)) = r 3 exp(-/3cr) = $^ 



(48) 



is just the Juttner distribution. 

For Bose-Einstein equilibrium distribution we have to = — 7 2 — j 2 R where R 
is defined in eq.(33). It can be seen from eq.(45) that in general 



for any equilibrium distribution 

5 Outlook 

There is the unique generalization of the relativistic diffusion equation (3) to an 
arbitrary metric [17] describing a diffusion in the presence of gravity (some other 
generalizations are discussed in [37] [38]). This is the diffusion which preserves 
the mass shell g^ u PnPv = m 2 c 2 (g^ v denotes the metric tensor). In [22] we 
have shown that an addition of a friction is necessary if the particle energy is to 
be bounded in time (then it equilibrates). The diffusion equation (3) (without 
friction) is explicitly Lorentz invariant as follows from cq.(6). However, the 
equilibrium measure cannot be Lorentz invariant if it is to be normalizable (in 
a finite volume) and if the mass shell condition (1) is to be satisfied (there is 
only one invariant measure on the mass shell, this is dpp$ 1 ). The equilibration 
to the Juttner or Bose-Einstein distribution takes place in a preferred Lorentz 
frame [39] [40]. In a general frame described by the four- velocity w v we should 
write the equilibrium distribution in the form 



where (by convention) in the rest frame w — (1, 0, 0, 0) and = A 1 (a 4 ^ — f° r 
k = 1,2,3). In general, all the parameters (/3,fj,,w) can depend on the position 
x (see [41] for such a framework in the description of the quark-gluon plasma). 

There is the unique way to generalize the Kompaneetz equation for a photon 
distribution function N(x,p) to an arbitrary frame on the pseudoriemannian 
manifold in such a way that in the rest frame on Minkowski space-time its 
linearized version coincides with the relativistic diffusion (9). The generalization 
of eq.(36) reads 



cxp(-2ft) = <&ei 



-l 



n E (w,n) = ( exp((3w u p v + (3w v ' ii„) - 1 



(49) 




(50) 



where 



D T = d T -Y = d T -ri a p a P °d j . 



(51) 
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Here, r are the Christoffel symbols and D is the covariant derivative on the 
cotangent bundle describing the geodesic motion. L flu (x) are the generators of 
the Lorentz group at the point x on the manifold and p are the coordinates of 
the cotangent bundle (the momenta, see [22]) 

Ljk = -i(g 3 iPkd l - g k iP 3 d l ) (52) 

and 

Loj = -igji\p\d l (53) 

The lhs of eq.(50) describes the well-known collisionless Boltzmann equation 
[42]. The operator A# (2) is the same as the one in the rclativistic diffusion 
equation (3). It does not depend on the metric (it is invariant under diffeomor- 
phisms of the momenta). In eq.(50) we suggest that the collision term resulting 
from the Compton scattering (which could be calculated on the curved mani- 
fold following the elementary derivation in [1]) has a trivial metric dependence. 
In the rest frame on the Minkowski space-time eq.(50) coincides with eq.(36). 
In an arbitrary frame N T has the Bose- Einstein distribution ue{w,^) as the 
equilibrium measure. 

In general, in astrophysics we have moving reference frames and non-zero 
gravitational fields. The effects of gravity seem to show no observational con- 
sequences unless the photons move in an electron gas around large compact 
massive objects. We have been concerned with the evolution at a large time 
in this paper. In astrophysical applications the behavior of the diffusion at 
large time (apart from its final effect: the equilibration) has not been discussed. 
However, with the increasing sophistication of the astrophysical observations the 
gravitational attraction as well as the large time effects of Compton scattering 
may give important information about the sources of the the CMBR distortions 
and the spectrum of light coming from compact stars. 
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